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Abstract 



We prove that the unique solution to the Yang-Yang equation arizing in the context 
of the thermodynamics of the so-called non-linear Schrodinger model admits a low- 
temperature expansion to all orders. Our approach provides a rigorous justification, 
for a certain class of non-linear integral equations, of the low-temperature asymptotic 
pi ^ expansion that were argued previously in various works related to the low-temperature 

behavior of integrable models. 



Introduction 

The non-linear Schrodinger model refers to a quantum field theory-based model for interacting bosons that is 
equivalent, in each sector with a fixed number of particles, to a one dimensional gas of bosons interacting through 
a 5-function potential with strength c. The so-called impenetrable Bose gas which corresponds to c = +00 was 
introduced and solved by Girardeau |7 1. Then Lieb and Liniger 1 14] introduced the model at general c and used the 
so-called coordinate Bethe Ansatz so as to construct eigenfunctions and eigenvalues associated with the underlying 
spectral problem. The completeness of their construction was however only proven much later by Dorlas IS. 

A method for studying the thermodynamics of the model was proposed in the seminal paper of Yang and Yang 
[16^ in 1969. These authors gave an integral representation for the free energy / of the non-linear Schrodinger 
model 



^ f, r, -£21. d^ 
/ = Jln[l+e .].- 



(0.1) 



The above representation is valid for any value < c < -1-00 of the coupling constant just as any value of the 
chemical potential h. The dependence on the coupling constant and the chemical potential is entirely contained in 
the function e which corresponds to the unique solution to the non-linear integral equation 



e{A)^ -h- ^ f K{A-^)ln[l+e-'-T^]-dfi with K{A) = _ _ . (0.2) 



^2 + c2 
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We shall refer to the above equation as the Yang- Yang equation. 

The method proposed by Yang and Yang underwent additional developments which allowed its implementa- 
tion to other integrable models such as the XXZ spin- 1/2 spin chain. For the latter model, the thermodynamics 
were characerized by a solution to a hierarchy of non-linear integral equations |6, 15 1. Later, in the nineties, there 
appeared an alternative method fTTl [T2l for characterizing the thermodynamics of spin chains. It allowed one 
to express the free energy of a spin chain in terms of the Trotter limit of the largest eigenvalue of the so-called 
quantum transfer matrix. In the case of the XYZ model, the latter eigenvalue was shown to be fully characterized 
in terms of a solution to a single non-linear integral equation |9| quite similar to the Yang- Yang equation (10.21 ). 
These further works on the thermodynamics of integrable models clearly point out the Yang- Yang equation as a 
prototypical form of a non-linear integral equation arizing in the study of thermodynamics of integrable models. 
Even if the resulting representations for the free energies, much in the spirit of (10.11) . may appear implicit, they 
still allow one to obtain the explicit behavior of thermodynamic quantities in the low-temperature regime. For this 
one first has to compute the low-J asymptotic behavior of the solution to the Yang- Yang equation (or its analogue 
for other models) and then built on it so as to extract the asymptotic behavior of the integral representation for 
the free energy such as, eg (10.11 ). Such calculations were considered by numerous authors and for many different 
models 012 13111113. 

The developments mentioned in the previous paragraph were not rigorous. For instance, the method set in 
'69 by Yang and Yang was formal and has been brough to a satisfactory level of rigor only 20 years later by 
Dorlas, Lewis and Pule [4] within the framework of large deviations. Likewise, the works IH HI |9l [TOl [131 solely 
proposed consistent algorithms allowing one to compute, order by order, the terms arizing in the low-J asymptotic 
expansion of the solution to non-linear integral equations of Yang- Yang type. However such algorithms always 
rely on the a priori hypothesis that the low-temperature asymptotic expansion exists in the first place and that the 
remainders remain stable under various operations. In this paper, we adress the question of building a rigorous 
approach to proving the existence of the low-temperature asymptotic expansion of the solution to (10.21) . The main 
result of this paper is: 

Theorem 0.1 Given a positive chemical potential h > and a positive coupling constant c > 0/or any integer 
^ e N, the unique solution to the Yang-Yang equation admits the low-T asymptotic expansion 

t 

£W = 2]£2iWr2*^ + 0(72^+2). (0.3) 

/(:=0 

Moreover, for any < a < c/2, the functions s^k ctre even and holomorphic in the open strip Sa where 

Sa = [z^C ■ |3(z)| < a\ . (0.4) 

Furthermore, the remainder in (10.31 ) holds uniformly in Sa and is stable in respect to a finite order X differen- 
tiation. 

The paper is organized as follows. In section [T] we review several known properties of the solution e(/l) to (10.21) 
and then establish some uniform in T bounds on the location of its zeroes. Then, several properties of integral 
operators that appear handy in the course of the proof of theorem lOHl are given in section |2l Finally, in section [3l 
we prove theorem lOTT] 

1 Several properties of the solution to the Yang- Yang equation 

We start this section by recalling, for the reader's convenience, Yang- Yang's proof for the existence of solutions 
to (10.21 ). We also provide a simple argument so as to justify the uniqueness of solutions. We then push further the 
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analysis of the objects involved in the proof of the existence of solutions what allows us to provide T-independent 
bounds for the zeroes of e(A). These bounds will play a role in the next section. 



Theorem 1.1 / 76/ The Yang-Yang equation admits a unique solution e{A) such that \_e{A) — /l^] € ^°°(R) fi (9(C). 
This solution is even, strictly increasing on ] ; +oo [ and satisfies to the bounds e'{A) > A for any A e ]0; +oo [. 
Lastly, for h > 0, s admits a unique zero 'q on ] ; +oo [. 

Proof — 

The Yang- Yang equation can be recast as the fixed point equation X[v](/l) = v{A) for the functional X : 

L°°(]R) ^ L~(R) given by 



X[/](^) 



-h - — ^ K{A-^^)\n[\+t-^]A^i. (1.1) 



Indeed, let v correspond to a fixed point of X, then it is readily seen that £{A) = v{A) + A? solves the Yang-Yang 
equation. The matter is that X has at most one fixed point. Indeed, given /, g e U°(K) 



X[/]W-X[g](^) 



-as— 



11/ - gllL"(R) 2) 



1 -Fexp|^max(||/|L,||g|U)) 



As a consequence, ||X[/] - XL^llli^^g-, < \\f - gW^o^^^^, so that the sharpness of the inequality ensures the unique- 
ness of a fixed point in L'"(]R). Note that two distinct solutions to the Yang-Yang equation would give rise to two 
distinct fixed points of Thus, there is at most one solution to the Yang Yang equation as well. 

We now build a sequence converging to the fixed point of Following Yang and Yang the sequence Vn{A) is 
defined as vq{A) - 0, v\{A) = -h and Vn+i{A) - £[v„]{A) for n > \. The sequence v„(/l) is such that 

• v„(/l) is a decreasing sequence in L°°(]R) of even functions of A; 

• v„(/l) is smooth and v',(/l) > for any A e]0; -i-oo [ . 

These properties are readily seen to hold for « = 1 Now, assume that these properties hold for some « € N*. Since 
Vn £ L'"(R), one has that £.[v„]{A) is also in L°°{R). Derivation under the integral theorems guarantee that X[v„](/l) 
is smooth. Since the integrand is negative, it is clear that V2 - vi < 0. Else, for « > 2 one has 

Vn+liA) - Vn{A) - Jk{A- //)| j — ^} • ^ < , 

and, using that v'„{A) is odd, one has for /I e ] ; -i-oo [ 



vL,(^) 



+ CO 

„l-i-er 1+er 



since, on the one hand v'„(jj.) > on [ ; -i-oo ] and on the other hand K{A - fi) - K{A + fi) > for {A,fi) € (R^f. 
As, A^ v„(/l) is increasing, it follows that v„(/l) > v„(0), so that 



v„+i(0) > L(v„(0)) with L(x) = -h - T J^(^)ln[l+e-'^-7]-^. 



(1.4) 
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The function 

L(x)-x = -h - T f *:Ou)ln[eT +e-^]- — . (1.5) 

R 

is strictly decreasing. The dominated convergence theorem ensures that 

lim L(x) - -h meaning that lim [L(x) - x] - -oo . (1.6) 

X—>+ca X— >+oo 

Moreover given any decreasing sequence x„ such that x„ ^ -oo, one has that the increasing sequence of functions 

- In l^e^ + e~^ j^r()u) is bounded from below by the integrable function - In [e^ + e~^ j^r(ju). Hence, by the 
monotone convergence theorem 

lim r-i.Ou)ln[e^ . e^^]^ = f = .oo . (1.7) 

R R 

Thus Wmx^+oo [L{x) - x\- +oo. Since L(x) - x is continuous on K. and strictly decreasing, it admits a unique zero 
-Zh on R. Moreover, since L(0) < -h and L{x) - x is strictly decreasing, it follows that z;, > for ^ > 0. One can 
also check that the function L is strictly increasing on R. Thus, 

lim L{x) - -h -h> L{x) for any x e K. . 

As a consequence, one has -h > L{-Zh) ~ -Zh- Then, a straightforward induction shows that v„(0) > -Zh- Since 
v„(/l) > v„(0) > -Zh, one has that v„(/l) is a decreasing sequence of measurable functions that is bounded from 
below. Thus v„(A) — > v(A) pointwise to a measurable function v(A). Taking the pointwise hmits in the relations 

Vn(A)>-Zh and Vn(A)-Vn(-A)^0 (1.8) 

ensures that v{A) satisfies to the same properties. As the sequence v„ is decreasing, = vo{A) > v{A), and thus 
V e L°°(R). 

It is then enough to apply the monotone convergence theorem so as to get that v(A) corresponds to the fixed 
point of X in L°°(R). Indeed, for any A, one has that 

li^hndi) = K{A-n)\n[\ +e —] (1.9) 

is an increasing sequence of positive measurable functions which, due to v„{A) > v„(0) > -Zh, is bounded from 
above by a L'(R) function 

2 

h„(M) < *:(0)hi [1 + e-'-^] . (1.10) 
Since hn converges pointwise, one has that 

h„(M)-^= K(A-ix)ln[l + e-^]-:f . (1.11) 

R R 

Taking the pointwise limit in v„+i(/l) = X[v„](/l) shows that v is indeed the fixed point of X- 

Finally, it follows from the integral representation given by v{A) = X[v](/l) that v is holomorphic in the strip 
Sc- By using this information and deforming the integration contour R in X[v](/l), it is readily seen that v is 



4 



holomorphic in S2c- Hence, by immediate induction, it is entire. A derivation under the integral sign ensures 
that v'(A) > for /i € ] ; +00 [. Hence, setting e(A) = + v{A) we get the sought unique solution to the Yang- 
Yang equation. Moreover, we immediately have that the latter satisfies e'iA) > on [ -17 ; 17], e(0) < -h and 
e(/^) > -Zh + A-^. Thus lim^_^+oo siA) = +00, so that e admits a unique zero ^on ] ; +00 [. Hence, by eveness, it 
has thus two zeroes ±^on R. 



Proposition 1.1 Let h > 0, then there exists To such that for any T e [0;Tq] the unique zero q of e as defined in 
theorem \L1\ satisfies V/z <'q <2 ^Jw, where w e ] ; +00 [ corresponds to the unique zero of 

, X , 2(X + C^) / 2c 

u(x) - X - h- arctani 1 + — yx . (1-12) 

n ^ c ' n 

Proof — Evaluating the Yang- Yang equation A - pleads to 

q^ = h + — \ K(,q- +&- — ]■ dn . (1.13) 

R 

The second term is positive meaning that 'q > ^/h, this uniformly in T € [ ; -1-00 [. 

It thus remains to obtain the upper bound. It was established in the course of the proof of theorem 11.11 that 
e(/l) > A^ - Zh where zh solves -Zh - L{-Zh) with L defined in (11.41 ). Thus "q < ^/zh and the upper bound will 
follow from an estimate on Zh- For this purpose, we recast the defining equation for Zh in the form 

^fzh +00 

Zh = h+ ^ (Zh - fi^)K(M) ■ ^ + VhiT) where Vh{T) = K(jj)ln[l + e-^] ■ dfi (1.14) 



It is readily seen that VhiT) < rin2. Explict integration of the second terms in (11.141 ) shows that Zh solves the 
equation tL>{Zh) = VhiT)- The function a> is smooth on ] ; -1-00 [ and strictly increasing 

co'ix) - 1 - -arctan(— ) > 0. (1.15) 

n ^ c ' 

Since a>{G) - -h and Oi){x) ~ 4c ^/x/n as x ^ +00, it follows that co admits a unique zero w on ] ; -1-00 [. Moreover 
a» is a C°°-differomorphism from ] ; -1-00 [ onto ] -h ; -1-00 [. Since VhiT) > 0, one has that Vh{T) e]-h; +00 [ so 
that Zh = (^~\VhiT)). Then, the Taylor expansion of a>~^(t) to the first order around t = and the estimates for 
VhiT) ensure that Zh = w + 0(T). In particular, as co~^ is strictly increasing, it follows that there exists Tq such 
that w < z/, < 4w on [ ; To ]. As a consequence, < 2 ^Jw uniformly in T € [ ; Tq ]. ■ 

2 Some properties of a linear integral equations driven by the Lieb kernel 
2.1 Invertibility of the Lieb operator 

Lemma 2.1 Let K{A) be the Lieb kernel defined in (10.21 ) and / € ([ —a ; a ]), then 



a a 



J f(jj)K{A-^i)d^idA\ < J \f(M)\dfi- jK{A)dA, (2.1) 
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Proof — 
One has 



a a a-ji 

Jf(p)K(A-^i)d^dA\ < Jd//{|/Ou)| J K{A)dAY (2.2) 



-a -a -a-fi 

The integral involving the kernel K can be decomposed as 

a-fi a -a a 

J K{A) -dA = J K{A) -dA + J K{A) ■ dA - J K{A) ■ dA . (2.3) 

-a-/! -a -a-/! a-fj 

Yet, for any ju € [-a; a], one has that 

-a a a+ii 

J K{A) -dA - J K{A) -dA = J [K{A) - K{A - fi)] ■ dA < . (2.4) 

—a—fi a—fi a 



This is evident for a > > 0, since then [K(A) - K{A -//)]< for any A> 0. The case when -a < ju <0 follows 
since the orientation of the integral produces a sign and [K{A) - K{A - fj.)] = [K{A) - K{A + |//|)] > for i > 0. 
As, the last two terms in (12.31 ) produce negative contributions, one gets 



J K(A)dA < jK(A)dA. (2.5) 



-a-fi 



Lemma 2.2 Let a e [0;+oo[ . The integral operator Kjln acting on L^iX—a ;aY) with the integral kernel 
K{A - jd) given in (10.21 ) is a compact trace class operator. The logarithm of the Fredholm determinant of I — K/ln 
is well defined and given by the convergent series 

K T laK{0) ^ 1 r A d"T 

I \ \ \ \ . ^2.6) 



lndet[_„;„][/--J- — - 2^-\\ \K{Ta-Ta,,)- 

n=l „ a=l 



{2nr 



where in each n-fold integral periodic boundary condition t„+i = t\ are being undercurrent. Finally, the operator 
I — K/2n is invertible and the kernel R'-'^XA, ji) of the resolvent operator R^"^ can be expressed in terms of the below 
series of multiple integrals 

+0O " (1-1 

R^"\A,ii) = K{A-ii) + y\ K{A-Ti)-r\ K{Ta - T,+i) • K{Tn - p.) ■ — ^ . (2.7) 



Given any < a < c/2, the series representation (12.71 ) is convergent in respect to the sup norm on Sa X Sa (with 
Sa as in (10.41) ). Moreover, for any (A,p) e S^, with < a < c/2 one has the estimates 

a 
-a 

and the resolvent is an analytic function of{A,p) on the strip Sc/2 X ■Sc/2 
Finally, the function {a,A,p) i-^ R^"\A,p) is on R^. 
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Proof — 

The estimates for any of the n-fold multiple integrals occuring in (I2.6I) - (I2.7I ) can be obtained by using the 
majoration (12. II ). We treat the case of the n-fold cyclic integral occuring in the series expansion for In det[/-/r/27r]: 



a a 

d'V r „, . „, . d"T 



Y\K{Ta-Ta^l)--^^ < KiO) J K{n - T2) . . . K{Tn-l - Tn) ■ 



{2nf 



- HI • / ^(ri-r,)...i^(r„^,-r„^i)-d"-ir < . (2.9) 

-a 

It is readily checked that 

r d/i 

J^W^ - 1 meaning that ||ii:|lLi([ 2;: . (2.10) 

R 

As a consequence, the series representation for In det [/ - K/2n] is convergent. The same types of estimates aslo 
ensure the uniform convergence in respect to the sup norm topology on compact subsets of M? of the series (12.71 ) 
defining the resolvent. 

Lastly, the fact that (a, i,//) ^ R^"\A,ij) e C\M?) can be readily inferred from the Fredholm series represen- 
tation for R'-''\A,ij.), det[_Q,;Q.] [/ - K/i2n)] and the fact that det[_Q,;Q,] [/ - K/i2n)] t for any given a e [0 ; +oo [. 



2.2 Positivity property of Lieb-like kernels 

Lemma 2.3 Let I — O be a trace class integral operator on L^([ -«;«]) with < a < +00 whose integral kernel 
is continuous and satisfies 

0{A,n) - 0(A,-fi) > and 0{A,^) ^ 0{-A,-fi) for any n,Ae]0;a[. (2.11) 
Assume that the series representation 

" n-l 

Ro{A,p) = OiA,p) + j 0{A,Ti)-Y\{0{Te,Ti^i)}-0{T„,fi)-d"T (2.12) 

for the resolvent kernel Ro{A,p) of the inverse operator I + Rq to I — O is convergent in L°°([-fl; a]^). Then, the 
resolvent kernel Ro{A, p) satisfies 

Ro{A,p)-Ro{A,-p) > for any e ]0;a[ . (2.13) 

Proof — 

We consider the series representation for Ro{A, p)-Ro{A, -p) and show that every term in this series is positive. 
Due to the hypothesis of the lemma regarding to the kernel 0{A, p), it is enough to show that 

n-l 



r 0{A,Ti) W 0{Te,TM) ■ [0{Tn,p) - 0{Tn, -p)] ■ d'V 
-a ^=1 

= j [0{A,Ti)-0{A,-Ti)]Y\[OiTe,TM)-0{Te,-TM)]-[0{Tn,p)-0{Tn,-p)]-d"T . (2.14) 
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For n - I, one has that 



a a a 

J 0{A,T)[0(T,n)-0{T,-^)}dT = J 0{A,T)[0{T,n)-0{T,-^)}dT + J 0{A,-T)[0{-T,fl)-0{-T,-^)]-dT (2.15) 
-a 

Above, we have spUt the integration contout, into the parts [ ; a ] and [ -a ; ] and have changed the variables 
T i-> -r' in the integral over the segment [-a;0]. It then remains to use the reflection invariance of the kernel 
0{-T, -fj.) = 0{t,iu) so as to get the representation for n = 1. Assume that (12.141 ) holds for some n. The ultimate 
integration over t„+i occuring in the n + 1-fold integration can be recast in much the same way as for « = 1: 

a a 
J 6>(t„,T„+i)[0(t„+i,//) - 0(T„+l,-yU)]-dT„+l = J [0{Tn,Tn+l) - 6>(t„, -T„+l)]-[6>(T„+i,yU)-C?(T„+l, -//)]-dT, 



Once this operation is perfomed, the integration in respect to ti , . . . , t„ is of the form taken care of by the induction 
hypothesis. The result is proven for « + 1. 
It then follows that 



Ro{A,n) - Roi^>-l^) = 0{A,^i) - 0{A,-n) 

n-\ 



2 r [0{A,tO - 0{A, -Ti)] ■ f] [0{Te,TM) - 0{t,,-tm)] ■ [0(t„,//) - 0(t„,-//)] • d'V . (2.16) 
Since every term in this series is positive, the claim follows. ■ 



2.3 Existence and uniqueness of the dressed energy 

In the Bethe Ansatz litterature, the dressed energy refers to the unique couple (eoiA \ q), q) solving the conditions 

S(){A I q) - f K{A - u)so(u \ q) ■ ^ = A^ - h where q is such that Eo{q \ q) - . (2.17) 

Lemma l2!2l ensures the existence and uniqueness of solutions to all integral equations of the type 

iI-K/27r)-g = f , forany f e L\[-a ■,a]) . (2.18) 

However, the unique solvability of (12.171) is a more involved issue in as much as the endpoint of integration q is 
also subject to a non-linear equation. In fact the statement about the existence of a solution to the system (12.171) 
has been made at several places in the litterature mEHHHllTOl altough, to the best of the author's knowledge, no 
proof has ever been given. In the proposition below we prove the unique solvability of the system (12.171 ). 

Proposition 2,1 Let eo{A \ a) denote the unique solution to the linear integral equation 



a 



sq{A I a) - \ K{A-n)so{n\a)-^ = A^-h. (2.19) 



Given any a, the function A ^ eq{A \ a) is analytic on the strip Sc/i- 

Also, there exists a unique ^ e ] ; +oo [ such that E(){q | ^) = 0. Moreover, there exists open neighborhoods 
Oq ofq and Vq ofO in M. such that the function a so{a \ a) is a -dijfeomorphism from Oq onto Vq. 
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Proof — 

Let R^"\A,ij) be the resolvent kernel associated with I -Kiln understood as acting on L^([ -a ; a ]). By lemma 
21 it is an analytic function of {A,ii) e Sc/2 x •Sc/2- Since 



£o(^ \a) = A^-h + J R^"\A,fi) .(i/-h)-^, (2.20) 



The statement about the analyticity of A eo(A \ a) on Sc/2 follows. Moreover, since (a,A,fi) i-> R^"\A,fj.) is 
and the integration in (12.201 ) goes through the compact interval [ -a ; a ], it follows by differentiation under the 

integral theorems that eoC'^ I ck) is continously differentiable in respect to a. 

In order to prove the existence and uniqueness of q, we study the variations of the function a i-> Eo{a \ a). For 

this, we first obtain a representation for its derivative: 



-^£o(a I a) = 4l^o(^ I a) U=q, + 
da oA aao 



— £0(a I a) = —£oi^ I a) \A=a + TT— eo(^ I ^o) \ao=a ■ (2.21) 



It is readily seen through integrations by parts that the last two derivatives satisfy to the integral equations 

[{l-^)-da£o]i^\a) = ^^^j^[K(A-a) + KiA + a)] (2.22) 
[(/ - £) • dA£o]{^ I «) = 2A- ^^^^[k{A -a) - K{A + a)] . (2.23) 

These are readily solved in terms of the resolvent kernel 

^so{^,a) = ^^^[R^<^\A,a) ^ R^'^\A,-a)] (2.24) 

a 

j-^soiA,a) = -f2|i^[/?(-Ua) - /?(-U-a)] + 2i + j Ri'^\A, f^)^df^ . (2.25) 

-a 

As a consequence, 

a 

—soia I a) = f2^^L^/?('»)(a, -a) + 2a + f [R^"\A,iu) - R^"\A, -//)]// • ^ (2.26) 
da TT J n 



It follows from lemma 1231 that the integrand of the integral appearing in the rhs is strictly positive on ] ; +oo [. It 

follows from (ITTl) that R^"\A,^) > 0. Thus —so{a | a) > as soon as eo{a | a) > 0. 

da 

Suppose that there exist aq > such that £o{q, q) - 0. It then follows from (12.261) that a i-^ e()(a | a) is strictly 

increasing in some open neighborhood ]q-r];q + M[ of q with 77 > 0, M > 0. Let M correspond to the largest 

real with this property and suppose that M + +00. One has that e()(a | a) is strictly increasing on ] ^ ; ^ + M [. 

Hence, by continuity E(){M + q\M + q) > Thus, one has that 

-^£o(a la) > , (2.27) 

da \a=q^M 

contradicting the definition of M. It thus follows that £0(0' I a) > on ] ^ ; +cx3 [. This reasoning shows that 
a £o(a I a) has at most one zero on [0 ; +00 ]. It thus remains to show that it has at least one. It follows from 
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(12.171 ) that eo(0 \ 0) = -h < 0. It thus remains to show that eo(a | a) > for or large enough. For this, recast (12.201) 
as 



£o(a \ a) = a'^ - h + J [R^''\a,fi) + R^"\a,-^i)](M^ - h)dn - J [R^''\a,^)+R^"\a,-n)]{h- A^)dn . (2.28) 

Clearly, the first two terms are strictly positive and behave at least as when a ^ oo. It follows from the 
bounds obtained in lemma fL2\ that the last term is, at most a 0(a). Thence, £o(a \ a) — > +oo. The function 

a 1-^ Eoia, a) being continuous, it admits at least one zero on ] ; +oo [ in virtue of the intermediate value theorem. 

Finally, since a i-^ e(a \ a) is smooth and [de(Qf | a)/da]a=q > 0, by the local inversion theorem for smooth 
functions, there exists open neighborhoods Oq of q and Vq of in R such that E{a | or) is a C°°-diffeomorphism 
from Oq onto Vq. ■ 



3 The low-temperature expansions 

In this section, building on the results of the previous two sections, we prove theorem 10.11 We start by proving 
the existence of the asymptotic expansion at order 0, what represents the hardest part of the theorem. Then, by 
using a convenient rewriting of Yang- Yang's equation (10.21 ). we show that the knowledge of the existence of an 
asymptotic expansion up to order 2^ + 2 as given in (10.31) guarantees the validity of this asymptotic expansion (10.31) 
to the order 2{ + 4. 



3.1 Asymptotic expansion to order 

Theorem 11.11 guarantees that the unique solution to the Yang-Yang equation is entire and has a unique zero ^ 
which, in virtue of proposition II. 1[ satisfies 'q > uniformly in T. Since s'iA) > 2A for A > 0, it follows that 
|£'(^)l > V/z uniformly in T. As a consequence, by the Bloch-Landau theorem, there exists an open neighborhood 

^^h\e'(q)\ h 

Uq c 2)^ ^ radms S > — > — such that e:Uq^ Dq^s (3.1) 

is a biholomorphism. Above, Da,r] stands for the open disk of radius 77 centered at a. In the following, we shall 
write for the associated local inverse of e, ie e~HO) - ^- This being established, we continue by decomposing 

? 

-T J i^(i-^)ln[l+e-^].| = J s(M)KiA-^)^ + T(°)(i) + ^(-)(i) . 

where 

^^-\A) ^ -T J^(i-;,)ln[l+e-^].^ and ^^^A) - -T J ^(^-//)ln[l+e-^].^ . (3.2) 

There, we have defined = e~\ ]-6 ■,6[), Ji^ = ] 5 ; -5 [ ) and we agree upon / ^ R \ U Ji^]. Note 
that both Jq and Jtq have a positive orientation. 
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Since e is even and strictly increasing on ] ; +00 [, it follows that \s{A)\ > 5 for any A e . As a consequence, 
for any /I e R, one has 

[^(■^^i)! < rin[l+e-7]x < rin[l+e-?]. (3.3) 

R 

Thus 'V^°°\A) only generates exponentially small corrections in T. The eveness of e allows us to bring the inte- 
gration in "V^"^ solely onto Then the change of variables ji = e^^{s) leads to 

T^'W ^ -T f('^<-'-''"<'>'"f-'"^"Wnin[l.e-¥].^. (3.4, 

J I e' o s~\s) I ^ ^ In 



Since 



inf \s'os-\s)\ > inf \£'{s)\ > inf |£'(5)| > Vh , (3.5) 



— ,1^ — 

it follows that, for any i e R, 



+00 

|^(0)(^)| < 4r2^^^ J ln[l+e-]d. = r^J'^"^ 



Hence, the Yang- Yang equation can be recast in the form 



(3.6) 



e(A) = A^ - h + jK{A-fi)e(M)-^ + ^^^^\A) + 'V^"'\A) (3.7) 
-? 

with the remainders V)(^) + V°°)(^) = 0(7^) being understood in the L°°(R) sense. By proposition O 
there exists constants w and To such that < 2 uniformly in T e [ ; To ] . It thus follows from lemma 12.21 
that / - Kiln understood as an operator on L^( [ -^;^] ) is invertible with inverse / -i- R'-'^^/Itt. Moreover, as 
is bounded for T € [0 ; Tq ], it follows from (12.81) that the resolvent kernel R^'^\A,fi) is bounded on R^ by a T 
independent constant. 

Thus recaUing that £o{A \ 'q) refers to the unique solution to equation ( 12.191 ) ( cf proposition 12.11) . it follows 

that 

£{A) - so{A\q) + r'^'\A) + r^^\A). (3.8) 
The last two function are given by 

? 



r(«)(i) ^ ^(")(i)_J^ /?<^\A,//)n^("V) with either u = or u = oo 



(3.9) 



As the resolvent kernel R^'^\A,^) is bounded on R^ by a 7 independent constant, it follows that 

V\-m = whereas = 0(e"T) = 0(7^) . (3.10) 
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This means that the zero q of e{X) solves the equation 

£® = eo(?l^ + 0{T^) = 0. (3.11) 

It has been show in proposition 12. 1 I that there exists a unique ^ € ] ; +oo [ such that so(q | ^) = 0. Moreover, there 
also exists open neighborhoods Oq of q and Vq of in R such that ni : a i-^ £o{a | a) is a C^-diffeomorphism 
from Oq onto Vq. Since for T small enough, one has 0{T^) € Vq, it follows that '5'= m-\-0{T^)) is well defined. 
Moreover, the Taylor expansion of around implies that '5" = q + 0{T^). 
Note that one has the decompositon 

? ? ? 

J £(ju)K{A - )") • ^ - J mK{A - • ^ + J e(M)[K{A - n) + K{A + y")] • ^ ■ (3-12) 

-? ~q 1 

We have already established that e{A) = £o{A \ q) + 0{T^) with the 0{T^) being uniform in i e R. Since eo(^ I ?) 
is uniformly bounded in some neighborhood containing q and "q, knowing that 'q = q + 0(T^) ensures that the 
second term in (13.121 ) is a again in the L°°{M.) sense. This means that, in fact, £{A) satisfies to 

£(A) = A^ - h + f K{A-n)e{n)^ + ^^^^\A) + 'y^°°\A) , (3.13) 
-<? 

where 

? 

'y(0)(A) ^ + J s{p)[K{A-ii) + K{A+ii)\-^. (3.14) 

1 

Again, •y^^/l) + 'V^'^XA) = 0{T^) in the L"(R) sense. Thus, by acting with the operator /2n on both sides 
of (13.131 ) and using, again, the estimates (12.81 ) we get that s{A) = so(A \ q) + 0(T^). The analytic properties of the 
resolvent kernel R^'^\A,fj.) obtained in proposition 12. 1 1 guarantee that, given any a € ]0; c/2 [, eo(A \ q) is indeed 
analytic on the strip Sa- As e is entire, it thus follows that the reminder has also to be holomorphic on Sa- Finally, 
it is readily seen from its explicit expression that the reminder is stable under finite order /l-differentiations. This 
last result can be as well derived from Cauchy's integral representation using that the remainder is holomorphic 
and uniform. 



3.2 Existence to all orders 

We now establish the existence of the asymptotic expansion to all orders in £ by induction. Thus assume that 
equation (10.31 ) given in theorem lOT] holds for some £. 

Observe that, in some vicinity of 0, the local inverse £~^{A) of e{A) admits the integral representation 

s-\A) = £ 4t^T^ fo'' ^ ^ ^o,<5/2 . (3.15) 
J £(z) -Aim 

For any dUq one has e{z) € dDo^s- Thus, for any A e Do,a/2 one has \£{z) - A\ > d/2 so that the integral repre- 
sentation (13.151) is indeed well defined. Moreover, this representation is well fit for providing a low-T asymptotic 
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expansion for the function in the neighborhood of as soon as the one for e is known. In particular, it implies 
that admits a low-T asymptotic expansion to the same order as e. 

So as to push the asymptotic expansion of e{A) two orders further, we first recast equation (13.131 ) in the form 

EiA) = SQ(A\q) + I^^\A) + I^^\A) + J<~\i), (3.16) 

where 

<5 - _j 

J(i)(i) - -T r/^^^i^-Mjinri+e-Tl.^ where R(A,i,) - R^^\A,iu) + R^^\A,-h) , (3.17) 

J I £' O E'\S) J L ^ 2k 



and 



I^'\A)^ js(M)-RiA,f^)-^ J<-)(i) . -rj^(^,;.)ln[l+e-^].|. (3.18) 



J' 



Moreover, below, we fix a particluar value of a, < a < c/2 defining the width of the strip First, we obtain 
the asymptotic expansion of I^^\A). It follows from the Taylor integral expansion 



i^,£-'{s)) V s'- y i R{A,e-\s)) ] 
' o E ^(s) ^ r] os'^ \ e' o e Mi) Jls=o 



R{A,s-\s)) _ ^ s' d' \ R{A,E-\s)) -\ 
s' 



with 



griA, s) = — (1 - m) • — — -J— du 

(r-l)!J ds'^y e' o £ \s) )\s=tu 





that one has the representation 



rn 4^T^'''^'^'> d^' (R(A,E-\s))] f., , i,,, ds f , ds 

(2r)! 5^'"^'^ I e' o e 1(5) J 1^=0 J 2n j 2n 

''-^ -i. -6 

T 

As the asymptotic expansion of e is uniform in respect to differentiations and that R is smooth with bounded 
derivatives on the strip Sa x <Sa, it follows that 

\g2t+2{A,s)\ < s'^^^Cc AeSa and se[-6;S] (3.19) 

where Cf is a T independent constant. Thence 



a +00 

T J g2{+2{A,s)ln[\+e-^'^]-ds < ICfT^^'''^ J u^^"^^ ln[l +&-"]■ du = OiT^^-"^) (3.20) 





Furthermore, up to exponentially small corrections in T, one can replace the integration over [ -6/T ; 6/T ] with 
one over R. The corrections will, again, be analytic, uniform in A e Sa and stable under finite order differentiation. 
Thence, using the integral representation for the Riemann ^-function 

+00 

(l_2-i-'-)^(r + 2) = — ^ r f\n(l+c-')dt , (3.21) 
r(r+l)J ' 
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we get 



J(i)(i) = - y 2r2(-i)(l - 2-'-''-)a2 + 2r)|^ f^,^'"''/;^H + 0(r2(^-^)) . (3.22) 



r=0 



By composition of the Taylor expansion of R{A,fi) in // belonging to a neighborhood of ^, and of the asymptotic 
expansions of £^^{s) and l/(e' o £'^{s)) into powers of T one readily gets that there exist functions /p\A), p = 
\,. . . ,{ analytic in the strip Sa such that 

Again, the remainder is stable in respect to finite order differentiations, uniform and analytic in Sa- 

Clearly, J(")(/1) - 0(r~), with a remainder that enjoys the usual properties. Thus, it remains to estimate 
I^^\A). Inserting the asymptotic expansion of e{A) up to 0(7^^^'*'^^) corrections leads to 

? ^ ? 

I'-^A) - j R{A,fi)-£o(M\q)^ + 2^'"/ m,n)-£2p(M)^ + 0{T^'^') (3.24) 

q P=l q 

We were able to replace the uniform 0(T'^^'*'^) appearing in the asymptotic expansion of e by a 0(7^'^"'"^) since it 
has been integrated over an interval of length \q-q\ = 0{T^). This 0(7^'^^'^) is differentiation stable, uniform 
and holomorphic in Sa- Since ^= £^'(0), ^admits a low-T asymptotic expansion to the same order as £ The 
regularity of the integrands then ensures that there exist analytic functions rf'\A), o^i Sa such that 

( ^ t 
J] r^/^ r R{A,n) . S2p(M)^ = J] T^''rf\A) + 0(7^'-') . (3.25) 

The remainder in (13.251 ) is also readily checked to have the desired properties. Moreover, we do stress that one 
has gained a supplementary precision 7^ on the remainder since on the one hand the sum starts with a and on 
the other hand, one integrates over an interval of length \q-q\ = O(r^). 

The first integral in (13.241) needs however a separate attention. Tayor expanding the integral to the order 
(q- qY^^ leads to 

J RiA,iu) ■ soiji \q)-^ = Yj ^jfiy ■ ^.{R(^,f^)soiM I q)]^ ^ + 0{(q- q)'^^) . (3.26) 



(r+1)! dfi'- 

q 



Note that the r = terms of this expansion is zero since £o(^ | ^) = 0. Thus, this expansion only contains terms 
(q - qY^^ with r > 1. Furthermore taking into accound that ^ - q = 0{T^) and that ^admits an asymptotic 
expansion to order 0(7^^''"^), it follows that (c^- qY'^^ with r > 1 admits an asymptotic expansion up to 0(7^^"""^), 
ie with a better precision than "q. Again, the regularity properties of the integrand then guarantee that there 
exist holomorphic functions r^j \A), . . . , rf^, (A) on Sa such that 



q 

I 



2 ■ ■ ■ e+v 

R(A,l^) .so(p\q)-^ = T'''rf\A) + 0{t^^'-'^) . (3.27) 

p=2 
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Adding up together all of the asymptotic expansions that we have obtained leads to 

p=i 

where we agree upon /^\A) = 0. Hence, as claimed, e{A) admits a low-J asymptotic expansion up to the order 
Q^j-2(f+2)^ with coefficients being analytic functions of A on Sa- Moreover, clearly, the remainder enjoys the 
claimed properties. H 



Conclusion 

We have proposed a method allowing one to prove the existence of the low-temperature expansion to all orders 
in T of the unique solution of the Yang-Yang equation. In the course of the proof, we have also established the 
unique solvability of the system (12.171 ) for the couple dressed energy-endpoint of the Fermi zone (eo(/l I q), q) 
characterizing the zero temperature behavior of the non-linear Schrodinger model. Our method mainly utilizes 
the positivity of Lieb's kernel so that it should provide one with analogous results for similar non-linear integral 
equations describing thermodynamics of other integrable models. Also, we expect that upon a few modifications, 
our techniques would allow one to approach rigorously the question of the large-volume behavior of various 
thermodynamic quantities arizing in integrable models, such as, for instance, the 1 /L corrections to the energies 
of the ground and low-lying excited states. 
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